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Abstract 

In automatic R-parity conserving supersymmetric (SUSY) SO(10) models, 
the simplest way to accomodate realistic fermion masses is to demand that 



the light Higgs doublets are linear combinations of the {10} and {126} grand 
unified Higgs representations. We study the realization of this mixed light 
Higgs property ( MLHP ) consistent with doublet-triplet splitting in a minimal 
R-conserving SUSY SO(10) model. We then discuss predictions for neutrino 
masses and mixings in this model as well as its implications for proton decay. 



1 Introduction 



Supersymmetric (SUSY) SO(10) model has a number of desirable features that 
make it an ideal candidate to describe physics beyond the standard model in a grand 
unified framework. They are: 

(a) unification of all fermions of each generation into a single spinor representation 
restoring quark-lepton symmetry to weak interactions; 

(b) a natural implementation of the see-saw mechanising for understanding small 
neutrino masses, which in the minimal version of the model are of the right order of 
magnitude to explain the solar neutrino puzzle via the MSW[[|] mechanism; 

(c) a simple mechanism for explaining the origin of matter starting with a zero 
baryon and lepton asymmetry of the universe for temperature above the grand 
unification scale 0]. 

In this paper, we wish to discuss a subclass of SUSY SO(10) models which 
have another highly desirable feature: automatic R-parity conservation which leads 
to natural conservation of baryon and lepton number symmetries prior to symme- 
try breaking. As is well-known, this property is not present in the SUSY standard 
model nor in the SUSY SU(5) model |J, where extra symmetries have to be imposed 
by hand to ensure R-parity invariance. On the other hand, in the SO(10) model if 
all Higgs representations are chosen to have congruence number zero (such as 45, 
54, 210, etc.) and two (such as 10, 120, 126, 126, etc.), the R-parity symmetry is 
automatic. Two possible minimal models with this property are given below: 

Model A : The Higgs particles belong to representations {210}, {126} © {126} 
and {10}@. The role of {210} is to break SO(10) down to SU(2) L x SU(2) R x 
SU(4) C ; that of {126} (denoted by A) is to break SU(2) R x SU(4) C down to U(l) y x 
SU(3)c while at the same time giving heavy Majorana mass to the right-handed 
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neutrino (i/r) to implement the see-saw mechanism for neutrino masses; the role 
of {126} (denoted by A) is to cancel the A contribution to the D-term so that 
supersymmetry is maintained down to the weak scale; the role of {10} (denoted 
by H) is of course to break the SU(2)^x U(l)y to U(l) e . m . and generate fermion 
masses. The model is minimal in the sense that omitting any one of these multiplets 
will leave extra undesirable local symmetries at low energies or break SUSY at GUT 
scale. 

Model B : The Higgs particles belong to the multiplets {45} © {54} (denoted 
by A and S respectively), and {126} © {126} © {10}. Apart from being more 
economical compared to model A, this model has another difference from model A. 
i.e., here SO(10) is broken down to SU(2) L x SU(2) /? x U(l) B _ L x SU(3) C before 
breaking to U(1) Y x SU(3) C - 

It is worth pointing out that in general superstring models, the {126} © 
{126} (or higher) SO(10) representations do not emerge below the Planck scale after 
compactification|7J. Therefore, several SO(10) models discussed in the literature M 
have only {45}, {16 H + 16 H }, and {10}, where {45} © {16 H } © {16 H } are used to 
break SO(10) down to the standard model. In these models, the vacuum expectation 
values(VEV's) of {16//} © {16//} lead to R-parity breaking terms at low energies. 
For instance, a term of the form \l/r a \E'\E'r a \l/ h /M p i ( where ^ denotes matter spinor 
and \I/// denotes Higgs spinor ) will lead to B-violating terms of the form < v c H > 
(u c d c d c + QLd c )/M p i, which can lead to catastrophic baryon number violationF]. 
One must invoke additional symmetries to prevent the R-non-invariant terms. It 
may very well be that such symmetries emerge from superstring compactification. 
But this remains to be demonstrated. Moreover, with the above minimal set, it 
is impossible to break the SO (10) symmetry down to the standard model without 
including large Planck scale induced non-renormalizable terms in the superpotential. 
1 The existence of the u c d c d c terms was pointed out to R. N. M. by A. Farragi 
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Let us now discuss the question of fermion masses in the R-parity conserving 
SO(10) models. It is well known |J that if only a complex 10 Higgs representation 
contributes to fermion masses, then one gets the undesirable mass relations between 
leptons and quarks e.g. md/m s = m e /m^, which is a factor of 10 off compared to 
observations. One way to correct it is to have the Higgs doublets contained in the 



126 contribute to the fermion masses through their direct dimension 3 couplings. 



However, if we naively added a separate 126, for this purpose, there would be an 
extra pair of Higgs doublets at low energies which is unacceptable from the point 
of view of gauge coupling unification^^]. Moreover, the usual problem of doublet- 
triplet splitting will also get worse. It was pointed out in Ref. that there is a 



simple solution to this problem: the same 126, which contributes to the breaking of 
B — L symmetry and to the see-saw mechanism, can have its weak doublets acquire 
an induced VEV without leaving any extra light doublet and without any fine tuning 



provided there is a coupling between the 10 and the 126 in the superpotential via 
other Higgs multiplets that only have VEV's of order of the GUT scale. (Such 
couplings can arise, for instance, if there is a 210 multiplet in the theory.) The 
induced VEV then will correct the bad charged fermion mass relations. 

This property of inducing VEV for (2,2,15) of A can also be stated in another 
way. When SO(10) breaks down to the standard model at Mjy, there are two light 
doublets (say <p u and 0^). To obtain realistic quark lepton mass relations, they 
must be a linear combination of the doublets in {10} and {126}. In the rest 
of the paper, we will call this "mixed light Higgs property" (MLHP). It is simple 
to see that models with MLHP do not have the property of Yukawa unification 



widely discussed in recent literature [ED] but strict Yukawa unification is anyway 



not realistic for the first and second generation. This property, MLHP does impose 
non-trivial constraints on model building. For instance, maintaining MLHP while 
implementing the doublet-triplet splitting (DTS) is rather nontrivial in general. 



In the Dimopoulos-Wilczek scheme [HI for DTS, an example was constructed |T3 



where this property emerges consistent with some discrete symmetries. In fact, it 



was suggested in Ref. [13| that any new DTS scheme must have this property In 
this paper, we will demonstrate a simple SO(10) model which has this good light 
Higgs property and study the consequences of this general class of SO(10) models 
with the additional minimality criterion. 

This paper has been organized as follows: in Sec. 2, we argue in favour of 
the necessity of the MLHP in R-conserving SO(10) models; in Sec. 3, we discuss 
the superpotential of the model B and introduce the doublet-triplet splitting mech- 
anism with MLHP; in Sec. 4, we study the quark and lepton mass matrices and 
discuss its implications for neutrino masses and mixings; in Sec. 5, we discuss the 
implications for proton decay; Sec. 6 is devoted to an R-conserving SO(10) model 
where the Higgsino mediated contributions to proton decay is naturally suppressed 
while maintaining MLHP; in Sec. 7, we summarize our results and conclude; in an 
appendix, we discuss the minimum of the Higgs potential consistent with the desired 
symmetry breaking pattern. 



2 Mixed versus Pure Light Higgs 

In this section, we explore to what extent, it is an absolute necessity to have 
mixed light Higgs doublets in an automatically R-parity conserving SO(10) model. 
i.e., could the light Higgs below GUT scale be purely light Higgs (PLH) arising solely 
from the complex {10} ? It is obvious that if non-renormalizable Planck induced 
terms are not included in the superpotential, then PLH scheme will not work since 
it will lead to the bad fermion mass relation md/m s = m e /m^ already mentioned. 
However, once the non-renormalizable Planck induced terms are included, the result 
is not obvious. To see what happens, let us assume first that the SO (10) is broken 
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down to the standard model by {210} Higgs, via VEV's for the components (1,1,1) 
© (1,1,15) © (1,3,15). The possible non-renormalizable terms involving the matter 
superfields are of the form: ^^§H, fyty$ 2 H, etc. It is then easy to verify that 
fermion mass matrices will have the general form: 

M u = (h + ti + 

M d = (h + ti -f-f')K d , 

Mi = (h + ti + 3/ + 3/')« d . (1) 

Here, /, /', and ti are of order y/8n Mu/M p i; f is symmetric coming from an effective 
126 operator; /' and ti are antisymmetric coming from effective 120 operators. One 
then has the relation 

Mi = 2rM u - M d . (2) 

Taking trace of both sides of this equation, we get r ~ m b /m t or zero. The right 
hand side of Eq.2 is completely determined by the known values of the quark masses 
and CKM angles. These values have to be extrapolated to the GUT scale in order 
to test the sumrule in eq.2. The details of this extrapolation procedure is described 
in Sec. 4 and is applied to the model B. Here we use those extrapolated values of the 
parameters to study the validity of Eq.2. We find that generically, the electron mass 
comes out a factor of five or so too large compared to the observed value with no 
free parameter left to adjust. The muon mass is also in disagreement with the known 
value by about a factor of 1.5. Therefore, we feel that it is highly unlikely that a 
pure light Higgs possibility in the sense defined in this paper would be realistic. 

Similar arguments apply if the GUT symmetry is broken by a combination of 
{45} © {54}. Thus, we feel that the mixed light Higgs doublet property provides a 
better chance to get a realistic fermion spectrum in a minimal R-conserving SO(10) 
model. In the next section, we give an example of an explicit model where MLHP 
is realized and proceed to discuss its implications in subsequent sections. 
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3 Mixed Light Higgs in Model B 

In this section, we present a simple minimal SO(10) model, where the light 
Higgs doublets have the desirable property (MLHP) described in the introduction. 
As already mentioned , the key step is to have a term in the superpotential that 
couples the 10 Higgs with the 126 Higg(s) via Higgs fields that have VEV's of order 
of the GUT scale. If we do not allow for Planck scale induced non-renormalizable 
terms, then we need a 210 Higgs to achieve this goal; but as shown infl3|], achieving 
this together with doublet-triplet splitting is not very simple, although we did man- 
age to construct an example which is technically natural. In this paper, we will take 
the point of view that one should include non-renormalizable Planck scale induced 
terms and we will keep only the lowest order Planck induced terms. As we will see 
this allows for the construction of a rather simple model with MLHP. 

As usual, we assign the fermions to the 16- dimensional spinor representation 
of SO(10). We denote them by ^ a (where a = 1,2,3 stands for generations) and 
we use the following minimal set of Higgs bosons needed for complete symmetry 
breaking: 

(i) {45} © {54} (denoted by A and S respectively) to break SO(10) down to the 
left-right symmetric group SU(2) L x SU(2) fl x U(l) B _ L x SU(3) C ; 

(ii) {126} © {126} to break the SU(2) L x U(l) B -£ symmetry down to U(l)y while 
keeping supersymmetry in fact down to the Mjy scale; 

(iii) A single {10} (denoted by H) to break the SU(2) L x U(l) y down to U(l) e . m .. 

The superpotential of the model is chosen to consist of the following parts: 

W = W f + W. + W p , (3) 

where 

W f = h ab V a y b H + f ab y a y b A, (4) 
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(fjL H + XS)HH + /j s S 2 + \ S S 3 + fi A A 2 + /i A AA 



+A A AAA + A S (SAA + SAA), 
^AA*H. 



M, 



P i 



(5) 
(6) 



As noted in Ref. ||, if we show that the light doublets in the model are 
linear combinations of the doublets in {10} and {126} multiplets, then Wf can 
accommodate a realistic charged fermion spectrum for all generations. This is, of 
course, intimately connected with the question of the doublet-triplet splitting. 

It is easy to see (see Appendix), using the above superpotential, that the 
vanishing of F-terms at the scale My and v r (which are equal to fit low energy LEP 
data||), is guaranteed by the following choice of VEV's for the Higgs fields S, A, 
A, and A: 



< S > 
<A> 

< A > (1,3,10) 



3 3 3 3 
diag{\, 1, 1, 1, 1, 1, --, --, --, --)Mu, 

ir 2 <8> diag(b, b, b, c, c), 
< A >(i, 3 ,10)= v R . 



(7) 



Using this, we find that the mass matrix for the fermionic doublets in H, A, and A 
can be written as: 



A d 
A d 



V 



A u 






A u 

yUA + AaC 



/'2 





/'i 




/ 



where we have fine tuned /i# — (3/2) XMu ~ M w (assumed to be zero in writing 
Eq. (H)). In Eq. (|8|), the rows and columns denote the down- and up- type doublets 

contained in H, A, and A respectively. The AA 2 H entries are induced by the Planck 
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scale corrections to W p leading to \i{ ~ / M pi ( i= 1,2 ) which is of order 

IQ- 1 M V to IQ-^Mu. Note that the (1,3,1) VEV in < A > makes the entry A u H d 
different from AdH u (i.e., /ii 7^ ^2). 

It is easy to see from Eq. (|8|), that the light Higgs doublets (denoted by 4> u and 
4>d) are given by 



cosaH u + sinaA u , 
cosjHa + sinjAd, 



(9) 



where tana = /ii/(//A — Aac) and tawy = ^/(/M + Aac). It is important to note 
that in general a 7^ 7. We find this to be a rather simple way to get the light 
doublets with the correct group theoretical property at low energies, for fermion 
masses. Furthermore, we expect a and 7 to be much smaller than one so that the 
departure from strict Yukawa unification |jll|is small. 

The triplet mass matrix is a four-by-four matrix as follows: 



H 
A 
A 



V 



H 

2/x H + \M V 


qi(-b 2 + a'c 2 ) 
q2bc 



A 

qi(b 2 + a'c 2 ) 
yUA - Aa/3'6 







A 




/iA + ^A/3'b 








y"A + Aa7'& ) 



(10) 



In Eq. (pT0[) , the first three rows denote the anti-quark-type Higgsinos contained 
in (1,1,6) of H, A, A, and the last row denotes the same in the (1,3, 10) of 
A respectively; similarly, the first three columns denote the quark-type Higgsinos 
contained in (1,1,6) of H, A, A, and the last column the same in (1,3, 10) of 
A. And the primed symbols represent non-zero group theoretical factors and qt are 
proportional to /ij. It is certain that all eigenvalues are of order Mjj. This solves 
the doublet-triplet splitting problem. 
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We want to point out that the specific form of the superpotential W$ + W p can 
be derived by requiring invariance under a Z 2 symmetry, under which S, H, and A 
are even and A, A are odd. This symmetry for instance forbids the AA 2 H term. 
Coming to the matter part of the superpotential, Wf, the \I/\I/A term is forbidden 
but there is an allowed Planck induced term f^^AA which essentially plays the 
same role as the second term in Wf. The effective coupling in the mass matrix is 
then given by / vr\/8tt ' jM v \. 

Before turning to a discussion of the fermion masses in the model, we wish to 
emphasize two points: first, in discussing the fermion masses in any R-conserving 
SO (10) model, one must first ensure that the light Higgs doublets have the correct 
MLHP property consistent with the doublet-triplet splitting. The absence of color 
singlet Higgs doublets in 45 representation makes our mechanism a good starting 
point for model building; secondly, this property of the mixed light Higgs doublets 
is not trivial to ensure while keeping all color triplet Higgsinos heavy. For instance, 
if instead of {45}, we used a {210} to break SO(10), the doublet mass matrix 
becomes a four-by-four matrix since it includes the (2,2,10) submultiplet of {210} 
and several elements in this matrix must be engineered to zero value to attain 



the MLHP goalfll3f . Of course, one could double the number of {126} © {126} 
multiplets such that one set contributed to the light Higgs doublets while the other 
breaks B-L symmetry and the two remain totally separate. However, in this case, 
one must worry about the possibility of unwanted pseudo-Goldstone supermultiplets 
which spoil gauge unification. 
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4 The fermion sector and predictions for the neu- 
trinos: 



Let us now turn to the fermion mass matrices in this model. It was shown 
that in Ref. ||, the fermion mass matrices in this model are characterized by 12 
parameters, which can all evaluated given the six quark masses, three charged lepton 
masses and three CKM angles. The neutrino masses and mixings are then completely 
predicted. We repeat this discussion with two differences from Ref. ||. First, we take 
the effect of the superpartners on the running of the gauge and Yukawa couplings. 
Second, we take the effect of the top quark Yukawa couplings on the running of the 
nj. We will follow Naculich|14| below. The low energy superpotential for 



masses 
the model is given by: 



W = h u Q4> u u c + h d Q4> d d c + h e L(p d e c + fjuj> u <f> d , (11) 

where h u , h d , and h e are three- by-three matrices expressible in terms of the SO (10) 
coupling matrices h and / in Eq. (£|) as follows: 

h u = hcosa + fsina, (12) 
hd = hcos r y + fsinj, (13) 
h e = hcos'-y — Sfsiwy. (14) 

As emphasized earlier a ^ 7 is required. Otherwise all CKM angles vanish. 
Moreover, in a strict supergravity framework, where supersymmetry breaking is 
implemented in the hidden sector, at fi = M p i, we have = m^ d = TOq. Their 
extrapolation down to the electroweak scale is governed predominantly |15j by h u ^3 



and . Correct symmetry breaking pattern (i.e., tan/5 > 1) also requires that 
a and 7 be different from each other. The electroweak symmetry is then broken 
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radiatively so that, 

< (fP u > = vsinfi, 

< 0° > = vcosft. (15) 

The mass matrices at GUT scale can then be written as: (// = M v ) 

Mu = (h + 7)v, 
Md = (hn+J^v, 

Mt = (hn-3jr 2 )v, (16) 

where 

h = hcosasin/3; f = fsinasinfl; r\ = cotfl; r 2 = — — cotf3. (17) 

cosa sma 

This is now in the same notation as in Ref. 0.1n order to evaluate the matrices h 
and /, T\ and r 2 , we use the following sum rule derived in Ref. ||, i.e., 

M l = ±^M u - 3 Jl±IlM d . (18 ) 
Ti -r\ r 2 - ri 



We then take Tr Mi, Tr M t , and Tr M { to obtain r 2 and J"i[16]. The light 
neutrino mass matrix is given by the see-saw formula|3| to beH 



where 



M v = —M v dM~mM v d, (19) 



M v u = — — — — — Af M —M d , (20) 

r 2 - n r 2 - r x 



M v m = -\ 



n — 1 — 
-M u .1/,/ 



^2 — n. r 2 — J~i 



(21) 



with 



vsinasinfi 
R = . 

12 



Note that this light neutrino mass matrix defined at vr needs to be extrapo- 
lated to the weak scale; but since z/r decouples below vr, there are only some over 
all anomalous dimensions of the effective light- Major ana neutrino mass|1l7|. This 
effect is small and we will ignore it. 

In order to carry out the numerical analysis, we choose the following values 
for the running masses: 

m u = 5.1MeV; m c = 1.27GeV; m t (m t ) = 166GeV; 

m d = (8.9 + ds)MeV; m s = (.175 + ss)GeV; m b = 4.25GeV; 
m e = .51MeV; m M = 105.6MeV; m T = 1.784GeV. (22) 



The m t (m t ) is obtained by taking the CDF[TI| mean value of 174GeV for the m t 
pole. The symbols ds and ss are left free to be fixed by the sum rule in Eq. QlSf) 
along with r\ and r-i- The CKM angles are parameterized in terms of s±2, S23, and 
S13, with S12 = —0.221, s 2 3 = 0.043, and S13 = 0.0045 as our choice corresponding 



to the mean values from experiments |19 



We then extrapolate all masses to the SUSY breaking scale[T4"]. The extrapo- 
lation factors are defined as r]i = m,i(mi or lGeV) /m^susy)- They are 

r] u = 2.17; rj c = 1.89; Vt = 1; % = 2.16; 
7] s = 2.16; 7] b = 1.47; i] e = ^ = r] T = 1.02. 

In order to extrapolate from Hsusy to Mu, we need to know tan (3. We follow 



NaculichJT^] and assume tan(3 < 40 so that effects of all Yukawa couplings except 
that of the top quark can be ignored. In this limit, the top Yukawa coupling effect 
is accounted for by the factor B t = 0.88647. The GUT scale values of the various 
masses (denoted with bars) are given by 



m u — m u rj u A u B t : 
m d = m d r] d A d Bt; 



m c = m c n c A u Bf; m t = m t f] t A u B^; 

m s = m s r] s A d B°; m b = m b r] b A d Bl ; 
13 



m e = m e r] e A e B^; m M = m^AjQ; m T = m T T] T A e B® , (23) 

where A-factors are the contributions of the gauge groups to the extrapolation and 
are numerically given by (choosing fisusY = 170 GeV) 

A* = 3.21; A d = 3A3; A e = 1.48. (24) 

Some of the mixing angles are also extrapolated and one has (for ij = 13 and 23 ) 

Sij = SijB-f. . (25) 

In order to predict neutrino masses, we need the value of the parameter R = 
vsinasinj3 /vr. The value of sina and sin/3 are arbitrary, whereas v = 246 GeV 
and vr is fixed by unification of the gauge couplings. In the absence of the heavy 
particle threshold corrections, one has Vr ~ My » 2 x 10 16 GeV[]Il|; but as has 



been noted for the case of non-SUSY SO(10) models[20], the threshold corrections 
can easily introduce uncertainty of a factor of 10 ±1 in My and vr. We will therefore 
assume that vr ~ 10 15 to 10 16 GeV in what follows. 

In our input, the signs of the fermion masses can be arbitray. There are many 
posibilities. Below we give the results for the choice of signs for masses that satisfy 
all the constraints of the model: 



Case I : All masses chosen positive. In this case, we find ds = —1.581, ss = 
-0.03533, n = 0.00952, and r 2 = 0.00479. The Eq. (|TJ) for this choice of {n,r 2 } 
leads to the values of {m e ,fn^,fn T } = {0.0003377,0.0695543,1.18177} to be com- 
pared with extrapolated values: { 0.0003378, 0.0695548, 0,1.18177}. The predictions 
for neutrino masses and mixing for this case are given by 

M v = #{-0.679348, 43.2421, 704.332}, 
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Vi 



\ 



0.999238 -0.0381721 0.00819604 
0.0379246 0.998875 0.0284804 
-0.00927398 -0.0281479 0.999561 



(26) 



A natural value for R w 10 -14 to 10~ 13 (since vsinasinfi ~ 10 2 GeV) depend- 
ing on whether vr is My or My/10, we get itl Vt « 7 x 10~ 3 eV to 7 x 10~ 2 , and 
m v k4x 10 _4 eV to 4 x 10~ 3 eV. The 6 eT mixing angle in this case is rather small; 
but sin 2 29 efl « 5 x 10~ 3 , which is of right order of magnitude to resolve the solar 
neutrino puzzle via the MSW mechanism [ pj]] . 



Case II : For the choice of all signs for masses to be negative, we get a consistent 
fit to all charged fermion masses for r% = 0.00637635, and r 2 = 0.0187798. The 
predictions for the neutrino sector in this case are 



M,. 



V, 



i?{-0.910023, -21.4474, -547.607}, 



/ 



0.995493 0.0921015 -0.0225942 
-0.0939442 0.990303 -0.102345 
0.012949 0.104006 0.994492 



(27) 



In this case, both 8 eT and 9 etl mixing angles are outside the MSW two neutrino 



solution given by Hata and Langacker||21||. Therefore, if the solar neutrino deficit 
situation continues to remain as it is now, this solution will be ruled out. 



Case III : We have found another fit to the sum rule in Eq. fli~8|), for the choice 
of masses, where m c , m^, and m s < whereas all the remaining masses are chosen 
positive. The values of r\ and r 2 are: r\ = 0.0100082 and r 2 = 0.072028. The 
predictions for neutrino masses and lepton mixing for this case are given by 

M v = #{0.0662987,4.60097,-2041.97}, 
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0.722321 0.684518 0.0984244 
0.69126 -0.718838 -0.0736991 
0.0203028 0.121271 -0.992412 



(28) 



Again, here all (mass difference) 2 are outside the range of the small angle as 
well as the large angle MSW solutions to the solar neutrino puzzle. Again, this 
solution can be tested by the solar neutrino data. 



Case IV : This case is obtained by changing the signs of masses in case III. This 
corresponds to the choices r\ = 0.0175744 and r 2 = 0.00699614. The predictions for 
neutrino masses and lepton mixing in this case are 



V, = 



i?{-0.460917, -30.2621, -663.422}, 



/ 



V 



0.99973 0.0226327 0.00523126 
-0.0227722 0.999339 0.0283497 
-0.00458617 -0.0284611 0.999584 



(29) 



Here, again the mixing angles are outside the range required by the MSW 
analysis of the present solar neutrino data. 

Case V : In this case, all masses are chosen positive except the electron mass and we 
find ds = -1.581, ss = -0.03833, n = 0.00955, and r 2 = 0.00500. The predictions 
for neutrino masses and mixing are given by 



M„ 



i?{-0.945754, 42.8837, 714.726}, 



0.999542 -0.029543 0.00655202 
0.0294036 0.999359 0.0204401 
-0.00715168 -0.020238 0.99977 



(30) 
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A natural value for R m 10 -14 to 10~ 13 (since vsinasinfi ~ 10 2 GeV) depend- 
ing on whether vr is Mu or Mjj/10, we get m VT ~ 7 x 10~ 3 eV to 7 x 10~ 2 , and 
m,y ~ 4 x 10~ 4 eV to 4 x 10~ 3 eV. The 8 efl and er mixing angles in this case are 
rather small for the MSW mechanism to work. 



Case VI : If we choose all masses to be negative except the electron mass, we get 
a fit consistent with all charged fermion masses for r\ = 0.00656483, and r2 = 
0.0187798. The predictions for the neutrino sector in this case are 



M,, 



V, 



#{-1.02584, -23.7002, -616.621}, 



/ 0.993278 0.112199 -0.0284463 ^ 



V 



-0.114706 0.98706 -0.112052 
0.015506 0.114562 0.993295 



/ 



(31) 



In this case, both 9 er and 6 efl mixing angles are outside the MSW two neutrino 
solution [51]. Therefore, if the solar neutrino deficit situation continues to remain 
as it is now, this solution will also be ruled out. 

If we set aside prejudices towards mixing angles coming from solar neutrinos, 
then our v il -v T mixing angles are in the interesting ranges to be testable in the 
next generation of proposed acceleration v [J -v T oscillation experiments. In Fig. 1, 
we compare our predictions with the domains of Am 2 and sin 2 28 UfiUT angles to be 
explored in the proposed CERN and Fermilab experiments. 



5 Proton Decay 

One of the key predictions of grand unified theories is the life-time of the proton 

and its decay mode. In non-SUSY GUT models, the dominant decay of the proton 

arises from the exchange of superheavy gauge bosons and the operators responsible 
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for this have dimension six. The primary decay mode is p — > e + ir°. On the other 
hand, in SUSY GUT models, in addition to the above dim.-six operators, there also 
exist dim-five operators, and in simple SUSY SU(5) or SUSY SO(10) models, the 
latter graphs dominate. The resulting dominant decay mode is p —>■ T7 fJ- K + , which 
can be used to distinguish between the SUSY GUT theories from non-SUSY GUT 
ones. 

Proton decay in SUSY SU(5) model has been extensively studied and 
it has been established that in this case, one requires the superheavy color-triplet 
Higgsino (H 3 ) mass M^ a ^> Mu in order to be consistent with the existing lower 
bounds on the r v ^ x+ P^fl ■ We will show below that in the SUSY SO(10) model 
dim.-five proton decay operators receive contributions from two diagrams: one in- 



volving {10} Higgs and the other involving {126} Higgs and these graphs could 
interfere destructively, thereby reducing the effective p-decay amplitude. This in 
turn can relax the constraints on color-triplet Higgsino masses. 

The color-triplet Higgsinos that mediate proton decay are part of a four-by-four 
matrix, given in Eq. ([11]). One can always find two four-by-four unitary matrices V 
and U, for the triplet mass matrix Eq. ([TUP , such that 

(V% k (M T ) kl U l:j = MA r 

Then, the dim.-five operators at the GUT scale, which are to be turned into baryon- 
number violating four-fermion interactions by gaugino- or Higgsino-dressing at the 
electroweak scale, have the following common factor: 

<^abcd - <^o 2^ 77 ' y 6Z ) 

where the early Latin indices are family ones and run over 1,2,3; Co and fo represent 
some over all factors (Clebsch-Gordan coefficients). 
From the definition Eq. fll6|) , we obtain 

T 1 r 2 M u - M d 

h = ■ , 

v -r 1 + r 2 
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j = l - ri M u + M d 

v —r 1 + r 2 

The first point to note is that since for all our solutions in Sec. 4 ri + r 2 rs 5 x 1CT 2 
to 5 x 10~ 3 , the proton decay amplitudes receive an extra enhancement factor ps 20 
to 200 (depending on the cases) compared to minimal SU(5). However, unlike the 
SU(5) case, we have two separate contributions to the proton decay amplitude and 
we could hope to invoke parameters for which there is a cancellation. To see if this 
is possible, we choose the special case where b 2 = a'c 2 . For convenience, we define 

h a b 

Vab = -7— ■ 

Jab 

We, then, get 

^abcd — ^O^abjcd 77 

i=l M i 

(Note that one of the triplet pairs completely decouples from the proton decay 
amplitude.) We have r] u ^ r] 2 2 ^ i]i2- We, thus, have three equations involving the 
three mixing angles that characterize the matrices Vu and Vn and we could therefore 
expect a solution for which the proton decay is suppressed. It must however be 
pointed out that this does require a fine tuning of parameters. 

6 Further Suppression of Proton Decay Ampli- 
tude 

In the previous section, we showed that the minimal R-parity conserving 
SO(10) model of Sec. 2 tends to predict higher strengths for the dim.-five proton 
decay operators compared to minimal SUSY SU(5); however, unlike the minimal 
SUSY SU(5) model, here there is the possibility of cancellation if one allows fine- 
tuning among the mixings and masses for the color-triplet Higgsinos. It is, however, 
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worth emphasizing that the predictions for the neutrino sector and realization of 
realistic quark-lepton masses are logically independent of the proton life-time pre- 
dictions. The question, therefore, arises as to whether it is possible to suppress the 
Higgsino-mediated proton decay amplitude while at the same time keeping the dis- 
cussion of light Higgs doublets with correct group theoretical properties to give the 
mass matrix structure of Sec. 4. A discussion of how to suppress such amplitudes 
in a different class of SO(10) models has recently given in Ref. [23]. We present 
a different example that has the addition property of "mixed light Higgs doublet" 
property (as defined in the introduction). 

We extend the model by addition of one more Higgs in the {10} representation 
(denoted by H2). Note that automatic R-conservation is maintained in the model. 
We choose the superpotential to be 

+X A SA 2 + /i A AA + X A AAA + ^AA 2 H 1 . (33) 

M p i 

We require all /i ss Mv and the fine-tuning condition /ii 2 + X s (—3/2)Mu ~ Mw- It 
is then clear that the light doublet mass matrix in Eq. @ is preserved. 

In Table I, we present a discrete Z 3 symmetry under which all but dimension 
two terms are invariant, providing a symmetry basis for this idea. Due to the pres- 
ence of dim.-two terms, this symmetry is softly broken, which in the supersymmetric 
context means that even after SUSY breaking terms are included no hard dim.-four 
term would be generated with infinite coefficient that break the symmetry (thereby 
helping to maintain our conclusion). This symmetry has the implication that it 
allows the following matter couplings that will eventually lead to realistic fermion 
masses and mixings if we assign ^ a ^b to be u under the symmetry. 

W' m = h^^Hx + ^ a V b A 2 A, (34) 



where M = M p i/\/8tt. The second term leads to both {126}, {120} as well as {10} 
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type couplings. The mass matrices are therefore less predictive than before. 



7 Summary and Conclusions 

In summary, we have studied the question of how to get realistic fermion masses 
in automatically R-parity conserving SUSY SO (10) models. We have argued that 
an economical way to do this is to have the light Higgs doublets contain a piece from 



the 126 Higgs multiplet that is responsible for the breaking of the B — L symmetry. 
Strict Yukawa unification does not hold in these models [TT]. But, of course, it is well 



known that models with strict Yukawa unification do not address the problem of 
second and first generation masses. We have given two examples of such models and 
studied their predictions for neutrino masses and mixings as well as proton decay. 
We find that all the models are testable by the solar neutrino experiments as well 
as the proposed accelerator experiments. In the proton decay sector, the situation 
is less predictive than in the case of the minimal SU(5) model due to the presence 
of two different contributions. We have also given a model in which the higgsino 
mediated proton decay can indeed be arbitrarily suppressed. 



Appendix 

In this appendix, we will discuss the symmetry breaking of the SO(10) model 
down to the standard model for the choice of Higgs multiplets S {54}, A {45}, 
A {126} © A {126}. The relevant part of the superpotential is given by W s in 
Eq. (||). To this we add the Lagrange multiplier term g TrS so that we can carry 
out the variation of all elements of the ten-by-ten symmetric matrix representing 
{54} © {1}. We will go to a basis where S is diagonal without loss of generality. 
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Let us now write down the constraints implied by all F-terms being zero. We will 
look for solutions with 

3 3 3 3 
<S> = diag(l,l,l,l,l,l,-~,-~,--,-~)Mu, 

< A > = ir 2 ® diag(b,b,b,c,c), 

< A > v c„c = < A >77 C ;7 C = v R , 

< H > = 0. (35) 
From the equation F s = 0, 

2/i s Mu + 3X s Mlj - \ A b 2 + g = 0, (36) 

2fi a (~Mu) + 3A S (^M^) - \ A c 2 + g = 0. (37) 

Demanding that TrS = determines g as follows: 

g= -^-KMl + iA A (36 2 + 2c 2 ) . (38) 

The other constraints are from vanishing of F A , F<\ and F% respectively: 

- 2^i A b - 2X A bM v + X A x d 2 = 0, (39) 
-2fi A c + 3X A cM v + X A d 2 = 0, (40) 
H^d + Aa(^o^ + c)d = 0, (41) 

where Xq is a Clebsch-Gordan (C-G) coefficient. 

Using the above constraints, we find the doublet Higgsino matrix with quantum 
numbers (2,1,1) or (2,-1,1) (under SU(2) L x U(l) y x SU(3) C ) to be that given in 
Eq. (|8[). As already mentioned, it contains only one pair of light doublets whose 
bosonic partners will be used to break the SU(2)^x U(l)y. For the color-triplet 
Higgsino matrix with quantum numbers (1,-2/3,3) or (1,2/3,3*) to be that given 
in Eq. (fL0|). As already mentioned, it does not have zero eigenvalues. 
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The Goldstone modes are contained in the following mixings: 
(2,1/3,3), (2,-1/3,3*); 5*2,2,6! ^2,2,6) ^2,2,i5> ^2,2,15 
(2,-5/3,3), (2,5/3,3*); S 2>2 , 6 , A 2)2 , 6 , 
(1,4/3,3); A 1)1)15 , A 1)3 ,io 
(1,-4/3,3*); Ai,i,i 5 , A 1AI5 
(1,2,1); A lj3 ,i, A 1Al5 
(1,-2,1); Ai, 3) i, A 1)3 ,io 
(1j0,1); Si,!,!, ^.1,1,15, ^1,3,1, A 13 jq, Ai, 3 ,i . 

Now, only two mixings are left. They are 
(1,0,8); S 112Q >, v4i,i,i5, 
(3,0,1); 5 3i3 ,i, A 3 ,i,i. 

We have exhausted all the submultiplets contained in A. The other submulti- 
plets of S, A, and A cannot be mixed. The remaining submultiplets in S are 
(1,-4/3,6), (1,4/3,6*); S lX2& , 
(3,2,1), (3,-2,1); 5 3 , 3 ,i. 
Their masses are given by 

\2(i a + 6\ a CiMu\, (42) 

where c, are C-G coefficients. We have checked that their masses are of order Mu, 
using Eqs. ©-(0), 

The remaining submultiplets are contained in A or A, and they are unmixed. 
Their masses have the following form: 

I A* a + A A (a^ + Vic) I, (43) 

where x% and t/i are C-G coefficients. To confirm that all the submultiplets are 
heavy, it is sufficient to show that Xj in Eq. fl4*3|) cannot be Xq in Eq. ( |41| ) or that 



Ui in Eq. (|3) cannot be 1. From simple group theoretical consideration, yi = 
for the submultiplets in Ai i, 6 , and A 31 jq. For the submultiplets which have 

23 



the quantum numbers (1,3,-2,1), (1,3,-2/3,3), and (1,3,2/3,6) (under SU(2) L x 
SU(2) fl x U(l) B _ L x SU(3) c ), contained in Ai )3j i , the ratios|5| of x { are 1:1/3:1/3. 
The submultiplet with the quantum numbers (1,-4,1) in the (1,3,-2,1) has the value 
Hi = —1. Considering tensor indices of A 2 ,2,i5 and the (1,0,1) in Ai^i 5 , we find that 
Xi = for the submultiplets in A 2j 2,i5- The same is true for the corresponding 
submultiplets in A. Thus it is guaranteed by Eq. ( J4"1"D that all submultiplets whose 
masses are given by Eq. ([|3p are superheavy. Also, we have checked that except the 
Goldstone modes all submultiplets involved in the above mixings have superheavy 
masses. 
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Table Caption 

Table I: Z 3 symmetry quantum numbers for various fields, (w = e l2w ^ 3 .) 



Table I 



Fields 


z 3 


S 


UJ 


A 


UJ 


#1 


UJ 2 


H 2 


1 


A 


UJ 2 


A 


1 



Figure Caption 



Fig. f : This figure shows the present limits on z/ M z/ r oscillation parameters (Am 2 and 
sin 2 29) and future possibilities on two proposed experiments CHORUS at CERN 
and P860A at Fermilab. The solid vertical lines are the predictions of the minimal 
SO (10) model described in this paper for the six allowed parameter ranges, denoted 
as cases I through VI in the text. 



28 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9406328v3 



